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ABSTRACT 

We determined the distances to nine pulsars by parallax measurements using 
the NRAO Very Long Baseline Array, doubling the number of pulsars with ac- 
curate distance measurements. Broadband phase modeling was used to calibrate 
the varying dispersive effects of the ionosphere and remove the resulting phase 
errors from the phase-referenced VLBI data. The resulting parallaxes have a typ- 
ical accuracy of 100 microarcseconds or better, yielding distances measurements 
as accurate as 2%. We also report new proper motion measurements of these 
pulsars, accurate to 400 /xas yr~^ or better. 



Subject headings: astrometry — techniques: interferometric — pulsars 
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1. Introduction 

Accurate pulsar distance measurements are useful for many reasons. Combined with 
dispersion measures, they probe the electron content along different lines of sight through the 
local Galaxy. Combined with proper motions, they constrain the pulsar velocity distribution 
and hence the symmetry of supernova explosions. Combined with radio, optical, x ray, and 
UV flux measurements, they determine the absolute luminosities and radiation efficiencies 
of pulsars. For a thermally emitting neutron star, they even promise the possibility of 
directly constraining the neutron star radius and hence the nuclear equation of state at 
high densities. There are many different ways of estimating pulsar distances, including 
association with supernova remnants or globular clusters at known distance, or the use of 
a model of the electron distribution in the Galaxy together with dispersion measurements. 
But the most fundamental method — indeed, the only model-independent method — remains 
the measurement of parallaxes. 

Pulsar parallax measurements are difficult. Even the closest pulsars are 100 pc or more 
distant; thus parallax leads to an apparent annual motion of 10 mas or (typically) much 
less. With their steep spectra, pulsars are usually observable only in the low-frequency radio 
band, where the ionosphere introduces substantial temporally- and geographically-variable 
shifts in the phase of the signal, making precision interferometry particularly challenging. 

Before the completion of the Very Long Baseline Array (VLBA), only six interferometric 
pulsar parallaxes were published (Table 1). It is important to understand the very great 
difficulties these observers faced. First, in many cases they were hmited to rather short 
basehnes (< 500 km) and decimeter and longer wavelengths (A > 10 cm). The ionosphere 
and the small array size limited the potential accuracy to about 1 milliarcsecond (mas), 
unless a very nearby calibrator source was used. Some serendipitous measurements were 
possible. For example, the measurement of the parallax of pulsar B1451— 68 was about ten 
times more precise than the ionosphere would ordinarily have allowed with the available 
array because of the fortuitous presence of a calibrator only 6' from the pulsar, within the 
same primary telescope beam. 

A second problem for early observers was the limited number of telescopes available 
in early, ad hoc VLBI networks. With only three or fewer stations, it is very difficult to 
fully resolve "phase-wrap" ambiguities to determine which of a family of possible positions 
on the sky is actually occupied by the source. As shown in Brisken et al. (2000), an 
incorrect choice of phase wraps could lead to parallax errors much larger than the formal 
measurement errors. (By the time the last measurement in Table 1 was made — the parallax 
of B2021-I-51 — a number of VLBA stations was available. Although the observations used a 
different set of telescopes at each observing epoch, and large baseline errors in the telescope 
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positions had still not been resolved, it is notable that this early result is within 2a of 
our new measurement.) It is probably not coincidental that every repeated pulsar parallax 
measurement has yielded a smaller parallax (and greater pulsar distance). We suspect that 
the expectation of larger results in early experiments probably contributed to incorrect lobe 
choices. With the ten-station VLBA this is no longer a problem. 

Although we will not further discuss non-radio measurements here, for completeness 
we note that the high-energy pulsar Geminga (J0633+1746) has a three-epoch parallax 
measured by HST of tt = 6 ± 2 mas (Caraveo et al., 1996), and the radio quiet neutron star 
J185635-3754 has an HST-derived parallax of 16.5 ± 2.3 mas (Walter, 2001). 
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Table 1. Early interferometric pulsar parallaxes 



Pulsar 


Stations 


Baseline 

(km) 


(MHz) 


epochs 


TT 

(mas) 


Reference 


B0823+26 


3 


5300 


1660 


4 


2.8 ±0.6 


1 


B0950+08 


3 


5300 


1660 


4 


7.9 ±0.8 


1 


B1451-68 


2 


275 


1600 


9 


2.2 ±0.3 


2 


B1929+10 


2 


127 


480 


~7 


21.5 ±8.0 


3 


B1929+10 


4 


35 


2695 


5 


< 4 


4 


B2021+51 


4-6 


~ 6000 


2218 


4 


0.95 ± 0.37 


5 



References. — (1) Gwinn ct al., 1986; (2) Bailes et al., 1990; (3) Salter et al., 
1979; (4) Backer & Sramek, 1982; (5) Campbell et al., 1996 
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2. Objectives and Target Selection 

Our successful measurement (Brisken et al. 2000) of the parallax of PSR B0950+08, 
using a new calibration technique to remove the disruptive effects of the ionosphere from 
the interferometry data, prompted us to expand our program of parallax measurements to 
a larger sample. We had two primary goals. One was the immediate desire to increase 
the number of pulsars with well-measured distances. But a second goal was to improve our 
understanding of the new calibration technique, to test its limitations, and to measure the 
level of residual systematic errors. For this reason, our target selection was not yet optimized 
solely to produce an unbiased pulsar sample for statistical work on pulsar distances. Instead, 
we chose nine pulsars at a variety of declinations, of varying flux densities, and with varying 
pulsar-calibrator angular separations. We also continued our observations of PSR B0950-I-08, 
to produce an extremely over-constrained data set for studies of our errors. 

An important constraint on target selection was the need for a high flux density at 
1400 MHz. Based on our experience with B0950+08, we expected that we would be able to 
successfully do an ionospheric correction for pulsars with flux densities above about 10 mJy 
at 1400 MHz. We also preferentially targeted pulsars with estimated distances (from the 
dispersion model) of less than about 2 kpc. For this experiment, we only considered candi- 
date pulsars whose nearest VLBA calibrator source was no more than 5° from the pulsar. 
Larger separations lead to larger differential ionospheric and tropospheric effects, and will 
be more difficult to calibrate. The 5° limit was set after experiments in which we attempted 
to calibrate PSR B0950+08 against the International Celestial Reference Frame calibrator 
source 1004-1-141, at a separation of 6.9°. At this separation the differential propagation 
effects were too large to correct using the new calibration technique. We also rejected pul- 
sars at declinations below —10°, where high elevation (> 20°) observations cannot be made 
simultaneously at all VLBA stations. The pulsars observed are hsted in Table 2, together 
with their dispersion measures and flux densities from the Princeton pulsar catalog (Taylor 
et al., 1993). 
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Table 2. Pulsars observed 



Pulsar 


DM 


'S'moo 


Gate gain''' 


Calibrator 


^sep 




(pc cm~^) 


(mJy) 


(mJy) 






B0329+54 


26.776 


203 


4.08 


J0302+5331 


4°37' 


B0809+74 


5.7513 


10 


3.27 


J0808+7315 


1°19' 


B0950+08 


2.9702 


84 


3.09 


J0946+1017 


2°52' 


B1133+16 


4.8471 


32 


5.98 


J1143+1834 


3°15' 


B1237+25t 


9.2755 


20 


5.13 


J1240+2405 


0°51' 










J1230+2518 


2°11' 


B1929+10 


3.176 


41 


4.03 


J1945+0952 


2°20' 


B2016+28 


14.176 


30 


4.99 


J2020+2826 


0°38' 


B2020+28 


24.62 


38 


4.69 


J2020+2826 


0°36' 


B2021+51 


22.580 


27 


4.90 


J2025+5028 


1°30' 


J2145-0750 


9.000 


10 


1.91 


J2142-0437 


3°18' 



^The improvement in signal-to-noise by not correlating during off-pulse. 
■1-B1237-I-25 was observed with two reference sources. 
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Tablc 3. Calibrator positions and flux densities' 



Name 


a(J2000)'' 


5(J2000) 


CC 

'-'2267 

(mJy) 


Cd 

'-'8337 

(mJy) 


J0302+5331 


03'^02"*22^7354 


53°31'46".534 


200 


192 


J0808+7315 


08''08™16".4918 


73°15'11".980 


220 


269 


J0946+1017 


09''46'"35^0693 


10°17WM26 


390 


248 


J1143+1834 


ll''43™26^0706 


18°34'38".375 


230 


113 


J1230+2518 


12''30™14^0912 


25°18'07".139 


150 


175 


J1240+2405 


^2^40^47^9870 


24°05'14".184 


160 


92 


J1945+0952 


19M5"'15^9229 


09°52'59".576 


250 


198 


J2020+2826 


20'^20"*45^8719 


28°26'59".205 


90 


18 


J2025+5028 


20'^25"'24^9726 


50°28'39".550 


280 


165 


J2142-0437 


2l''12™36\8999 


-or3ri3".5i8 


65 


33 



^AU values listed are from the VLBA Calibrator List (Beasley et 
al, 2002). 

^Source positions are accurate to 15 mas. Better positions for some 
of these calibrators are now known. 

^Correlated flux density on 8000 km basehne at 2267 MHz. 

'^Peak source brightness of VLBA image made at 8337 MHz. 
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3. 



Observations 



Five epochs of observations for each pulsar were planned over the course of one year, 
enough to redundantly determine the proper motion and parallax even with a failed epoch. 
To maximize our sensitivity to the angular signature of parallax, the measurements were 
scheduled at times when the parallactic displacement was at a maximum or minimum. Ob- 
servations hence occurred at roughly three month intervals, although the exact dates varied 
with telescope scheduhng. Pulsar J2145— 0750 was not detected in either of its first two 
epochs. This non-detection is most likely because of the extreme, long-term scintillation 
modulation that affects the signal from this pulsar. It was dropped from the program, and 
the time allotted to other pulsars was increased slightly. About ninety minutes of on-source 
integration was obtained for each pulsar at each epoch. 

In order to make relative position measurements, each pulsar observation was interleaved 
with observations of a VLB A calibrator source with position known to better than 15 mas. 
The antennas nodded back and forth between the pulsar and calibrator with a cycle time 
of about five minutes. This interval is long enough to detect calibrator fringes yet short 
enough to allow unambiguous phase connection between calibrator observations. In the case 
of B1237-I-25 two calibrators were observed in a cycle, integrating for 100 seconds on the 
pulsar and 60 seconds on each calibrator. 

All observations were made with the ten-station NRAO VLB A. The 20 cm band (1400 
to 1740 MHz) was chosen as a compromise between falling flux densities and increasing 
resolution as frequency increases. This band also offers the wide fractional band widths 
needed for the ionosphere cahbration. The 338 MHz wide band available from the receivers 
was Nyquist sampled in eight 8 MHz spectral windows beginning at 1404, 1414, 1434, 1474, 
1634, 1694, 1724 and 1734 MHz. Pulsar gating was employed to increase the signal-to-noise 
of the pulsar measurements by disabling correlation during off-pulse. The gate gain listed in 
Table 2 is the expected signal-to-noise increase factor based on the pulse profiles measured 
at 1400 MHz. The timing data needed to construct the pulse arrival time ephemerides was 
collected by Andrew Lyne at Jodrell Bank. 



Data reduction was performed with the AIPS software package. Phase-referencing 
(Shapiro et al. 1979) was used to relate the pulsar's position to that of its calibrator. 



4. 



Data Reduction 



4.1. 



Calibration 
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Corrections for the structure of the cahbrator source were made. The absolute position shift 
caused by structure is less than 1 mas for most VLBA calibrators. Epoch-to-epoch position 
shift caused by differences in the baseline coverage is a smaller effect but is more important 
for relative astrometry. Phase-referencing cancels out most of the pulsar's phase errors such 
as those due to clock and baseline offsets and bulk propagation effects; however, propaga- 
tion effects that are path-dependent are not exactly cancelled and dominate the remaining 
phase errors. At 20 cm, the largest of these effects is due to the ionosphere. A technique 
for removing the effect of the ionosphere without independent calibration information was 
described by Brisken et al. (2000). It is briefly described here. 

The phase shift at frequency v due to a free electron column density of o^- is 

« (1) 
The wet troposphere also causes phase shifts, but with a non-dispersive dependence on i/, 

A0^°P° oc an.o ^, (2) 

where cTHaO is the column density of water in the atmosphere. Only the differential phase 
shifts at two antennas enter into the visibility phases. After phase-referencing, only the 
difference in visibility phase shifts due to propagation between the pulsar and calibrator 
remain. While the absolute phase shifts may be several thousand degrees, only the double- 
differences of their values are relevant to the position measurements. For the troposphere, 
we flnd that this amount is typically less than 45° of phase for calibrator-target separations 
of up to 5° at 20 cm. In contrast, the differential ionosphere can be up to 1200° on long 
basehnes. 

In addition to the unwanted phase shifts due to propagation effects, the geometric 
phase associated with the position of the pulsar (the quantity that is to be determined) also 
contributes to the pulsar's visibility phase. The point-like nature of pulsars simplifies the 
geometric phase to 

0^*^°™= -(£M + m^;), (3) 
c 

where (£, m) are the pulsar's offset from the correlation center in radians, c is the speed of 
light, and (m, v) are the components of the projected baseline perpendicular to the direction 
to the pulsar. The phase delays due to geometry and the troposphere are both independent 
of frequency (non-dispersive), making their disentanglement impossible for a given visibility. 
Visibility measurements spanning a large range of elevations can be used to separate out the 
troposphere, but this was not done due to the limited elevation ranges that we used. For- 
tunately, it is possible to measure the differential ionosphere strength with multi-frequency 
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data since the ionosphere's phase delay is frequency dependent (dispersive). The visibihty 
phase on a single baseline can be expressed as 

^Aiy + ^ + 2Tm. (4) 

where A incorporates all of the non-dispersive components (the geometry and the tropo- 
sphere), B is the strength of the differential ionosphere, and the integer n is the number of 
additional phase wraps. Once B is known, the distorting effects due to the ionosphere can 
be removed from the visibility. 

Station-based ionosphere solutions were determined by fitting station-based phases to 
Eqn 4. The station-based phases at each frequency band were determined by fitting the 
measured phases to a point source model representing the pulsar with solution intervals 
shorter than the amount of time that the ionospheric phase shift changes by about 30° 
(about 20 or 30 seconds at 20 cm). Solution intervals of 10 seconds or less were used on the 
brightest pulsars. Since all phases are relative, one station was chosen to be the reference 
antenna and was assigned zero phase. 

Unknown phase wraps complicate the fitting for A and B. This is because certain com- 
binations of the A and B parameters closely mimic an additional undetectable phase wrap 
due to the limited spanned bandwidth as is shown in Figure 1. The fitting was performed 
first on the Southwestern Antennas (Fort Davis, TX, Kitt Peak, AZ, Los Alamos, NM, and 
Pie Town, NM) where baselines are shorter and thus the geometric phases are smaller and 
the differential propagation effects are also less. An image was made without any iono- 
sphere calibration to determine the pulsar's position to better than 10 mas. The pulsar's 
approximate position was used to reduce the unknown non-dispersive phase to less than 90°. 
Since the tropospheric phase is usually less than 45°, the number of phase shifts, n, can 
be determined and the ionosphere can be calibrated away for the Southwestern antennas. 
The phase shifts associated with the fitted B parameters were removed from the visibility 
data. An image made with the ionosphere calibrated Southwestern antennas yielded pulsar 
positions to better than 1 mas, allowing ionosphere solutions to be found at and applied to 
the remaining stations. 
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onosphere Solution Ambiguity 




0.5 1 1.5 2 2.5 

V [GHz] 

Fig. 1. — The source of phase-wrap ambiguity. For certain ionosphere solutions, A* = 
2.0 rad/GHz and B* = 4.9 rad-GHz, the resultant phase across the observed band is almost 
indistinguishable from 2tt. Since 2n can be added to the visibility phases without observable 
consequences, an entire family of ionosphere solutions is consistent with the data. The eight 
8 MHz bands observed are shown as short vertical bars. The RMS deviation from 27r at 
these frequencies is only 0.6°. 
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Ionosphere Solutions on Short Baselines 
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Fig. 2. — The relative phase due to the different residual ionospheric content above two 
antennas that is left uncalibrated by phase-referencing. Kitt Peak (KP) was used as the 
reference antenna in this observation. The three shortest basehnes to Kitt Peak are shown 
above. Distances are: Fort Davis (FD) to KP = 744 km, KP to Los Alamos (LA) = 652 km, 
and KP to Pie Town (PT) = 417 km. These solutions are from the phase-referencing of 
PSR B1929+10 via calibrator J1945+0952, a separation of 2°20'. The gap in the data 
between 100 and 115 minutes is due to the observation of the bandpass calibrator. The larger 
gap on the KP-PT baseline is due to telescope failure at PT. The date of this observation 
was 1999 November 16. 
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lonosphere Solutions on Internnediate Baselines 
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Fig. 3. — Same as Figure 2 but on tlie three intermediate length basehnes to Kitt Peak. 
Distances are: Brewster (BR) to KP = 1913 km, KP to North Liberty (NL) = 2075 km, and 
KP to Owens Valley (OV) = 845 km. The bias in each frame is due to systematic differences 
in the line-of-sight electron density at the two stations. 
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onosphere Solutions on Long Baselines 
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Fig. 4. — Same as Figures 2 and 3, but on the three longest basehnes to Kitt Peak. Distances 
are: Hancock (HN) to KP = 3623 km, KP to Mauna Kea (MK) = 4466 km, and KP to 
Saint Croix (SC) = 4839 km. For these longer baselines the systematic bias is greater. 
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Differential ionosphere strengths, converted to their effect on visibihty phases at 1400 MHz, 
are shown as a function of time on all baselines to Kitt Peak, the reference antenna for this 
observation, in Figures 2, 3, and 4. The phase differences on the four shortest baselines 
shown (KP to FD, LA, OV, or PT) all have RMS values less than 1 radian, and are centered 
close to zero phase. The longer baselines often show large systematic offsets from zero phase. 
See for example the BR-KP and KP-MK baselines. Without correction a displacement of 
the pulsar from its true sky location results. The data in these figures is from the 1999 
November 16, observations of PSR B1929+10. The sun was 68° from the pulsar during this 
observation. Sunset occurred at Saint Croix during the observation and is probably respon- 
sible for the rapid change in ionospheric phase delay near the end of the observation on the 
KP-SC baseline. Observations made when the pulsar and sun are on opposite sides of the 
Earth show considerably smaller ionospheric phase delays but still exhibit fluctuations on 
the 10 minute timescale. B1929-I-10 is separated from its calibrator by 2°20'. 



Final pulsar images were made once the ionosphere calibration was completed. A sep- 
arate image was made for each 8 MHz band. The position of the pulsar was measured by 
fitting a Gaussian ellipse to the pulsar image. Data from a single 8 MHz spectral window 
produced images of the pulsar with signal-to-noise of about 40. The dominant source of 
the image noise was phase scatter due to the troposphere, making the final signal-to-noise 
ratio almost independent of the pulsar's strength. In most cases the position agreement 
between images made at two different frequencies was better than 0.1 mas. A misalignment 
of more than this amount prompted a detailed inspection of the calibration process and 
usually resulted in stricter elevation limits or the removal of a frequency band with severe 
radio frequency interference. Once image alignment was achieved, a final image was made 
incorporating the data from all useful frequency bands. The final position measurement was 
made from this image. The synthesized beam sizes for the final images were about 10 mas 
north-south and 4 mas east-west but varied slightly depending on source declination. 



We model the path of the pulsar in terms of its position at epoch 2000.0 (q;2ooo, (^2000), 
proper motion {/lon/J's), and parallax (tt) by 



4.2. 



Position determination 



4.3. Calculation of Parallctx and Proper Motion 




(5) 
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where t is time since year 2000.0. The functions fa and fs are the parallax displacements in 
right ascension and declination, respectively, for an object at a distance of 1 pc as a function 
of the time of year. The five parameters are determined by minimizing the function 
associated with the above model using an iterative linear inversion. Uncertainties in the fit 
parameters were based on the single-epoch position uncertainties. 

5. Results 

The proper motion and parallax results for the nine pulsars with new parallaxes are 
shown in Table 4. The distances to the pulsars derived from their parallax measurement 
are shown in comparison with their dispersion measure (DM) derived distances in Table 5. 
The uncertainties in DM distance arc obtained by assigning the DM a Gaussian probability 
distribution function (PDF) centered on the measured DM with a 40% variance. The DM 
PDF provides a way to account for the uncertainty in the dispersion model; it is not meant 
to reflect uncertainty in the the dispersion measures as the DMs of the pulsars studied in this 
sample are measured to better than 1%. Since the dispersion model is not smooth, a non- 
Gaussian distance PDF is obtained when the dispersion model is applied to the Gaussian DM 
PDF. The most probable distance is quoted along with errors bracketing the most compact 
68% likelihood region. B 1237+25 has a very high upper limit because of is to its position 
high above the galactic plane {z ~ 850 pc). The dispersion model provides poor constraints 
to its distance since very little ionized material is thought to be at or above this pulsar's 
height. Also listed in Table 5 are the pulsars' transverse velocities {v±) and the average 
electron densities (ng) along lines of sight. 

The derived parameters of six other pulsars with accurate parallax distances are shown 
in Table 6. Chatterjee et al. (2001) measured the parallax of B0919-I-06 using the VLBA. 
The proximity of an in-beam calibrator to the pulsar (6'sep < 20'), and their simultaneous 
observation, reduce the amount of uncalibrated ionosphere and troposphere. Further, GPS 
data was used to model and remove most of the remaining differential ionosphere. Single 
epoch uncertainties as low as 0.2 mas were attained in this manner. The other flve pulsars 
are millisecond pulsars with distances determined through pulse arrival timing. B1534+12 
has a timing parallax upper limit of vr < 1.8 mas. A more accurate distance of d = 1.08±0.15 
was inferred for this pulsar by monitoring the decrease in its binary period, Pfy. By assuming 
that General Relativity accurately predicts the orbital decay of this close neutron star binary 
via emission of gravitational radiation, the other contributions to Pi, from the motion of the 
system restricts its possible distance. (Stairs et al. 1998 & 1999). 

The distances to the pulsars in Tables 5 and 6 are plotted in Figure 14. The dispersion 
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measure distances to these pulsars are plotted against the parallax distances in Figure 15. 

The dispersion measure distances tend to be smaller than the observed parallax distances. 
This may be due to the local low electron density bubble not being included in the dispersion 
measure model. For more detail on the local bubble, see Toscano et al. (1999). 
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Table 4. Position, Proper Motion and Parallax Results 



Pulsar 


•^2000 




^2000 




(mas/yr) 




(mas/yr) 




TT 

(mas) 


B0329+54 


03''32'"59^ 


.3862 


54°34'43", 


.5051 


17.00 ±0, 


.27 


-9.48 ±0.37 


0. 


,94 ±0.11 


B0809+74 


08'*14™59", 


.5412 


74°29'05", 


.3671 


24.02 ±0, 


.09 


-43.96 ± 0.35 


2, 


.31 ±0.04 


B0950+08 


09''53™09", 


.3071 


07°55'36", 


.1475 


-2.09 ±0, 


.08 


29.46 ±0.07 


3 


.82 ±0.07 


B1133+16 


11^36™03", 


.1829 


15°5r09", 


.7257 


-73.95 ±0, 


.38 


368.05 ±0.28 


2, 


.80 ±0.16 


B1237+25* 


12^39"^40", 


.3589 


24°53'50", 


.0193 


-106.82 ±0, 


.17 


49.92 ±0.18 


1, 


.16 ±0.08 


B1929+10 


19''32"^13^ 


.9496 


10°59'32", 


.4178 


94.82 ± 0, 


.26 


43.04 ±0.15 


3 


.02 ± 0.09 


B2016+28 


20'^18"'03^ 


.8332 


28°39'54", 


.1527 


-2.64 ±0, 


.21 


-6.17 ±0.38 


1 


.03 ±0.10 


B2020+28 


20'^22"^37^ 


.0718 


28°54'23", 


.0300 


-4.38 ± 0, 


.53 


-23.59 ± 0.26 


0, 


.37 ±0.12 


B2021+51 


20''22"^49^ 


.8655 


51°54'50", 


.3881 


-5.23 ±0, 


.17 


11.54 ±0.28 


0, 


.50 ±0.07 



pulsar's epoch 2000.0 position, (q;2ooo, (^2000) is based on the position of the reference source 
used (see Table 3) and is accurate to about 15 mas. 

^B1237±25 data is based on its closest cahbrator, J1240±2405. 
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Table 5. Derived parameters 



Pulsar 


di)M 


4 


v± 






(kpc) 


(kpc) 


(km s-i) 


(cm ^) 


B0329+54 


1 4+0-3 

-■-•^-0.5 




951- 


0.0253 ±0.0030 


B0809+74 


U-<J-L--0.13 


Q .00+0.008 


1021^ 


0.0133 ±0.0002 


B0950+08 


0.16lg;g? 


n 262+°-°°^ 

U.ZUZ_o.o05 


36.6l°? 


0.0113 ±0.0002 


B1133+16 




0.35+°:°^ 


631+1 


0.0136 ±0.0008 


B1237+25 




0.851^^ 




0.0108 ±0.0007 


B1929+10 


A 1 7+0.06 


n qqi +0.010 
U.OOl_Q 


1631^ 


0.0096 ± 0.0003 


B2016+28 


n 7+0.3 
-0.3 


0.951^^ 


30^^ 


0.0146 ±0.0013 


B2020+28 


1 q+0.5 
^^•'-'-0.5 


9 q+1.0 
^•'-'-0.6 


256+^i^ 


0.0092 ±0.0030 


B2021+51 


1 9+0.5 
-■-•^-0.5 


1 q+0.3 

-'-•^-0.2 


1151}^ 


0.0113 ±0.0016 



''"Gaussian parallax uncertainties imply non-Gaussian distance 
uncertainties. The most compact 68% uncertainty interval is 
shown above. 

•I-The mean electron density to the pulsar, rie, is DMd~^. 
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Table 6. Derived parameters of other pulsars with accurate distances 



Pulsar 






v± 


Tie 


Reference 




(kpc) 


(kpc) 


(km s"-*-) 


(cm~^) 





109+i; 0.0148 ±0.0021 1 

671}^ 0.2319 ±0.0477 2 

4841^^ 0.0226 ± 0.0036 3 

ISltjj 0.0104 ±0.0014 4 

231^ 0.0146 ± 0.0040 5 

2Sl:f 0.0144 ± 0.0048 6 

35tl 0.0088 ± 0.0006 7 



J0437-4715 O.Utom 0.17;[5;[]^ 

B0833-45 O-Qltliy 0-28^0^5 

B0919±06 61^^ l-15+o;?6 

B1534±12 0.71^^6° l-OSlaM 

B1855±09 0.7toi 0.791°;?? 

J1713±0747 O.Stoi 0.9toi 

J1744-1134 O.nt^om 0.351H1 



References. — (1) Sandhu et al., 1997; (2) Caraveo et al., 2001; (3) Chat- 
terjcc ct al., 2001; (4) Stairs ct al., 1999; (5) Kaspi et al., 1994; (6) Camilo et 
al., 1994; (7) Toscano et al., 1999 
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Motion of B0329 + 54 




15 10 5 -5 

Aa [mas] 

Fig. 5. — The modeled path of B0329+54 based on parameter estimates (see text) and its 
four measured positions. The eUipses represent la position uncertainties. 
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Fig. 6. — The modeled path of B0809+74 and its five measured positions. 
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Motlon of B0950+08 
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Fig. 7. — Eight observations of B0950+08 made over seven epochs and its modeled path. 
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Motion of B1 133+16 
■2-10 1 2 



3 



3 - 



2 



CO 
D 



- 



2 



-3 - 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




2000.698 



J I I I I I I I I I I L 



- -3 

J I I I I I I I I I I I I I I I I I I : 



- 



-3 -2 



-1 1 

Aa [mas] 



3 



Fig. 8. — A plot showing the measured positions of B1133+16 at 5 epochs and its modeled 
path, both with the estimated proper motion removed. This pulsar's proper motion of 
375 mas yr~^ ranks as one of the highest. The proper motion was subtracted to allow the 
significance of the parallax to be visible. 



-25- 




Fig. 9. — A plot showing the measured positions of pulsar B1237+25 at 5 epochs and its 
modeled path, both with the estimated proper motion removed. The linear motion of this 
pulsar moves it 118 mas each year. 
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Fig. 10. — A plot showing the measured positions of B1929+10 at 5 epochs and its modeled 
path, both with the estimated proper motion removed. The proper motion was subtracted 
before plotting since its linear motion over one year amounts to 104 mas. 
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Fig. 11. — The modeled path of B2016+28 and its five measured positions. 
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Motion of B2020+28 in RA 
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Fig. 12. — The large distance to B2020+28 makes its parallax measurement the smallest 
of those observed. Its parallax fit is shown above with its proper motion subtracted. Be- 
cause this pulsar's single-epoch position measurements are comparable in magnitude to the 
estimated parallax, the parallax fit is better visualized one dimension at a time. 
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Fig. 13. — The modeled path of B2021+51 and its five position measurements. 
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t= 90 



B2021+51 n,=0-.0ri3 



B0329+54 n,=0.0253^9L.' 

B0809 + 74 n^^-0. 01 33 O 
■■ B1929.M-10 n =0.003 




B2020+28 n =0:0092 



""■0-B2O16+28 =0.0146 

VbI 855+09 n,=.0.0146 

^ B1534+12 n3=0.0104 

O ,,.-Gr^^1 71 3+0747 n,=0.0144 
J1744-:1 134 n,=0.0q88 



4= 



; B0950+08 n3=0.0113 
B1 133+16 n:,=0.0136 



B0919+O6 n =0.0226 



O B1 



J0437-4715;' n^=0.0148 
B0833-45 n,=0.2319 

237 + 25- n^=0. 01 08 

1 kpc 



2 kpc 



3 kpc 



Fig. 14. — A projection of pulsars with accurate distance measurements onto the galactic 
plane. Each pulsar is labeled with its name and estimated line of sight electron density, rig. 
The radial line through each pulsar is its distance error bar, representing its most compact 
68% confidence interval. 
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Fig. 15. — Dispersion measure derived distances compared against parallax distances. The 
solid error bars are for the nine pulsars studied in this work (listed in Table 5); the dotted 
error bars are for those pulsars in Table 6. Error bars in each dimension represent the 
most compact 68% confidence interval. The predominance of pulsars to the lower right of 
the dashed diagonal line suggest that the dispersion model systematically underestimates 
distances to nearby pulsars. 



-32- 



6. Uncertainties and Limitations 

Several error estimators were used to set the final parallax errors. Pulsar images made 
at each of the 8 MHz bands were fit for position. In almost all cases scatter in position 
was about 0.1 mas or less. For the brightest pulsars this scatter in position dominated the 
final uncertainty. The other single-epoch consistency check was making images at several 
different short time ranges (~ 30 minutes) and observing the apparent motion of the pulsar 
over the ~ 2.5 hours of its observation. The dominant motion was always in right ascension, 
never amounting to more than 0.3 mas. Excising low elevation data (< 20°) usually reduced 
this disagreement to less than 0.1 mas. The source of this wander is almost certainly differ- 
ential refraction in the troposphere and its effect scales with calibrator-target separation at 
approximately 0.05 mas per degree of separation, although it is very likely that the scaling 
is non-linear for separations larger than those probed with these observations. Based on the 
tropospheric limitations, calibrator separations greater than the 5° limit that was initially 
set would probably not generate single-epoch position measurements with accuracies better 
than about 0.3 mas. While station weather data was used to estimate tropospheric refrac- 
tion, more sophisticated tropospheric modeling techniques exist and should be employed in 
future observations. Correlator model deficiencies, such as incorrect Earth orientation or 
station positions or calibrator positions, can also cause relative motion between the target 
and calibrator. The success of geodesy esperiments at the 10 fias level using the same cor- 
relator model as was used for this work suggest that these effects are not significant for 
these observations. The post-fit residuals provide the final check on the uncertainties. The 
least squares fit quality parameter, x^/(2A^epochs — 5), was less than 1.5 for all nine pulsars, 
indicating that the single epoch error estimates are reasonable. 

Pulsar B1237-I-25 offered a unique test of the final uncertainties. The pulsar and two 
calibrators were observed together in a slightly longer three-pointing nodding cycle. This 
allowed the relative proper motions and parallaxes of each pair to be measured. Table 7 
shows the proper motions and parallaxes measured for each pair, which are in excellent 
agreement. 

The small calibrator-to-pulsar separation is not the only limitation to the ionosphere 
calibration. The fitting of the phases imposes sensitivity requirements. Phase solutions 
accurate to ~ 10° are required on timescales of 30 seconds or less for each frequency band. 

B0950+08 now has seven epochs of VLBA astrometric data spanning just over 2.5 years. 
Figure 7 shows these measurements and its estimated proper motion and parallax. The RMS 
fit residual in right ascension was 0.22 mas, and that in declination was 0.15 mas, similar to 
the uncertainty estimated by its single-epoch position uncertainties, typically 0.2 mas. Its 
distance, 262 ± 5 parsecs, is the most accurate pulsar distance measurement made to date. 
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Table 7. Calibration Consistency Test 



Calibrator 


Target 






TT 






(mas yr~^) 


(mas yr~^) 


(mas) 


J1240+2405 


B1237+25 


-106.82 ±0.17 


49.92 ±0.18 


1.16 ±0.08 


J1230+2518 


B1237+25 


-106.53 ±0.24 


50.21 ±0.23 


1.08 ±0.11 


J1240+2405 


J1230+2518 


0.10 ±0.23 


0.09 ±0.19 


0.13 ±0.11 



7. Conclusions 



The ionosphere cahbration method used in this paper successfully removed the effect of 
the ionosphere from the VLBA data for all nine detected pulsars, allowing sub-milliarcsecond 
parallax measurements to be made. The requirements for this technique to work can be easily 
summarized - station-based phases must be measurable to ~ 10° at each frequency band 
with 30 second or shorter solution intervals. For our observations, the limiting flux density 
at 1400 MHz is about 8 mJy. Observing with a larger bandwidth, dividing the bandwidth 
into fewer spectral windows, or adding larger VLBI stations to the array would allow this 
technique to operate on weaker pulsars. Up to 100 more pulsar parallaxes are within reach 
of the VLBA. 

The pulsar timing data used for gating the correlator was very kindly provided by 
Andrew Lyne of Jodrell Bank Observatory. The National Radio Astronomy Observatory 
is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. A National Science Foundation graduate student fellowship 
funded Walter Brisken. Additional funding was provided by Princeton University and by 
the NSF (AST-0098343). 
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